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We present the first correlated-electron theory of metam- 
agnetism in strongly anisotropic antiferromagnets. Quantum- 
Monte-Carlo techniques are used to calculate the field vs. 
temperature phase diagram of the infinite-dimensional Hub- 
bard model with easy axis. A metamagnetic transition sce- 
nario with 1. order and 2. order phase transitions is found. 
The apparent similarities to the phase diagram of FeBr2 and 
to mean-field results for the Ising model with competing in- 
teractions are discussed. 



In 1939 Becquerel and van den Handel reported 
the observation of an anomalous magnetization behav- 
ior in the mineral mesitite (carbonate of Fe and Mg) at 
low temperatures which they could not explain in terms 
of existing theories of magnetism. Therefore they sug- 
gested for it the name "metamagnetism" 0. Qualita- 
tively similar, but even more drastic magnetization ef- 
fects were later observed in many other systems of which 
FeCki and Dy 3 Al 5 0i2 (DAG) are well-studied prototypes 
||. These materials are strongly anisotropic antiferro- 
magnets where the spins are constrained to lie along 
an easy axis e. Under the influence of a strong mag- 
netic field along e they undergo a 1. order ("metam- 
agnetic"), rather than a spin-flop, transition from the 
antiferromagnetic (AF) to a paramagnetic state. To- 
day the term "metamagnetic transition" is used in a 
much wider sense Q, namely whenever x(H) nas a max- 
imum at some value H c , with M(H) being strongly en- 
hanced for H > H c S. Metamagnetism is then found 
to be a rather common phenomenon which occurs also in 
strongly exchange-enhanced paramagnets (e.g. B TiBe2, 
YC02), heavy fermion and intermediate valence systems 
(e.g. CeRu2Si2, UPt3), and Q the parent compound of 
high-T c superconductivity, La2Cu04. 

In this paper we will be concerned only with metamag- 
netism in strongly anisotropic antiferromagnets, where a 
spin-flop transition does not occur. Apart from the well- 
known insulating materials, e.g. FeCl2, there are also 
some good electrical conductors, especially UAi-^B^ 
(where A — P, As; B = S, Se), in this group ||. Theo- 
retical investigations of metamagnetism in these systems 



were sofar restricted to the insulating systems, several 
of which are known to have a very interesting H — T 
phase diagram (H: internal magnetic field, T: tempera- 
ture). In particular, it includes a tricritical point (TCP) 
at which the 1. order phase transition becomes 2. order 
(Fig. la) - a feature very similar to that found in 3 He- 
4 He mixtures HJ^]. Theoretical investigations of the mul- 
ticritical behavior are usually based on the Ising model 
with competing interactions Q, e.g. 



H — J S ' SiSj — J' S ' SiSj — H S ' Si 
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For J, J' > one has an AF coupling between the Z 
nearest-neighbor (NN) spins and a ferromagnetic cou- 
pling between the Z' next-nearest- neighbors (NNN). It 
was pointed out by Kincaid and Cohen Q that in mean- 
field theory a TCP as in Fig. la exists only for R = 
Z'J'/(ZJ) > 3/5, while for R < 3/5 this point sepa- 
rates into a critical endpoint (CE) and a bicritical end- 
point (BCE) (see Fig. lb). The latter behavior, espe- 
cially the finite angle between the two transition lines at 
CE and the pronounced maximum at the 2. order line, 
is qualitatively very similar to that observed in FeBr2 
||. However, the 1. order line between CE and BCE 
has so far not been observed - neither experimentally, 
nor even theoretically when evaluating (|l|) beyond mean- 
field theory |o|Q . Most recently, by measurement of the 
excess magnetization and anomalous susceptibility loss, 
de Azevedo et al. jl2| identified a strip-shaped regime of 
strong non-critical fluctuations in FeBr2 which is at least 
reminiscent of the critical line CE <-> BCE. This unusual 
behavior was then also found theoretically by Selke and 
Dasgupta |l3| who evaluated ([!]) in the limit of weak 
ferromagnetic coupling (J' /J = 0.2) and large Z using 
Monte-Carlo techniques. 

It is the purpose of this paper to go beyond an effective, 
classical Ising-spin model and to investigate the origin of 
metamagnetism in strongly anisotropic antiferromagnets 
from an even more microscopic point of view. Our ques- 
tion is: which is the minimal electronic correlation model 
for metamagnetism in these systems? Since the Hubbard 
model at half filling (n = 1) is the generic microscopic 
model for antiferromagnetism (both of itinerant and lo- 
calized nature) in correlated electron systems, we will 
study this model. For nearest-neighbor hopping of elec- 
trons on a bipartite lattice in the presence of a magnetic 
field it has the form 
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NN,a i 

- ^(ji + vH)n l<7 , (2) 

where operators carry a hat. In this paper we consider 
a half-filled band in which case fi = U/2 by particle- 
hole symmetry. From the exact, analytic solution in di- 
mensions d = 1 the (paramagnetic) ground state of this 
model is known to exhibit metamagnetic behavior, i.e. 
dx/dH > up to saturation [fL4| . The only other dimen- 
sion in which the dynamics of the correlation problem can 
be treated exactly in the thermodynamic limit is d = oo 
jL5j. In this limit, with the scaling t = t*/V~Z in (g), 
one obtains a dynamical single-site problem Jl6| which is 
equivalent to an Anderson impurity model complemented 
by a self-consistency condition fPij] , and is thus amenable 
to numerical investigations within a finite-temperature 
Monte-Carlo approach |Q. This approach has recently 
provided valuable insight into the physics of strongly cor- 
related electron systems, e.g. the Mott-Hubbard transi- 
tion M| an( i transport properties [M . The effect of the 
magnetic field H in (Q) for n = 1 was also studied |H|-p3[ . 
In particular, Laloux et al. j23] thoroughly investigated 
the magnetization behavior of the paramagnetic phase, 
assuming the AF order to be suppressed. For U = 3^/2t* 
they find a 1. order metamagnetic transition between the 
strongly correlated metal and the Mott-insulator at a 
critical field H ~ 0.2t*. Giesekus and Brandt [^lj also 
took into account the AF order. They considered the case 
where the field H orients the staggered magnetization 
m st perpendicular to itself. For isotropic antiferromag- 
nets this is indeed the arrangement with lowest energy. 
There is then no metamagnetism. 

To investigate how metamagnetism may arise in corre- 
lated electron systems with strongly anisotropic AF order 
we will work in the AF phase of the Hubbard model, too, 
but will constrain m st to lie parallel to H . In this way 
the existence of an easy axis e along which H is directed, 
such that e || m s t \\ H, is incorporated in a natural way. 
We will first discuss (||) in the limit of strong and weak 
coupling since this will already show that the appearance 
of a tri- or multicritical point is a delicate matter. 

1. Strong coupling: In the limit U ^> t the Hubbard 
model at n — 1 is equivalent to an effective Heisenberg 
spin-model, H Heis = JJ2nn J2 a =x, v ,z - # E 4 §*, 

with an AF exchange coupling J = 4t 2 /U. For this model 
Weiss-mean-ficld theory becomes exact in d = oo yielding 
the same results as for the Ising model (|l|). In this case 
the transition line in the H — T phase diagram has indeed 
the form shown in Fig. Ia - but it is of 2. order for all 
T > 0, i.e. Tt = 0. The behavior changes if we include 
the next term in the effective spin-Hamiltonian, of order 
t 4 /C/ 3 , involving two- and four-spin terms EJ]. It leads 
to an effective ferromagnetic spin coupling J'. Then we 
obtain precisely the Ising model (0) in mean-field theory 



with J = At 2 /U and J' cx {t/U) 2 J < J, and thus expect 
to find the transition scenario in Fig. Ib. However, since 
the expansion is in t/U <C I it is not clear down to what 
values of U/t this behavior can actually be observed. 

2. Weak coupling: For U -C t we expect the Hartrec- 
Fock approximation, a static mean-field theory, to give 
an, at least, qualitatively correct answer. Within this 
approximation we find a metamagnetic phase transition, 
too - but it is of 1. order for all T and U, corresponding 
to Fig. Ia, but with Tt = Tn- Hence the Hartree-Fock 
solution cannot describe the experimental situation. Ap- 
parently the location of the (tri-) critical point in the H - 
T phase diagram, and even the transition scenario itself, 
depends sensitively on the value of the electronic on-site 
interaction U . To study this point in greater detail we 
have to go to intermediate coupling. 

3. Intermediate coupling: In this interaction range we 
solve (||) numerically in d — oo. In contrast to Hartree- 
Fock this limit provides a dynamic mean-field theory. 
The local self-energy Yf an and propagator G^„, where 
the subscript n denotes the Matsubara frequency u n — 
(2n + X)-kT and a S (A, B) = (+, -) is the sublattice in- 
dex, are determined by two sets of dynamically coupled, 
self-consistent equations for G and S 1 16 1?],|25| : 



G" 



de- 



N°(e) 



(3) 



Here z^ n = iuj n + jjb — and the thermal 

average in is defined as a functional integral 

over the Grassmann variables Tp,ij)*, with (0)A a = 
J V[ip]V[ip*]Oexp{A a {ip,ip*,i:,G})/Z a , in terms of the 
single-site action 



-U / dr^f(r)V» T (r)^(r)^(r) , (4) 
Jo 

with Z a as the partition function. Furthermore N°(e) 
is the density of states of the non-interacting electrons. 
As the results do not much depend on its precise form 
we choose N°(e) = [(2t*) 2 - e 2 ] 1 / 2 /(27rt* 2 ). From now 
on we will set 2t* = 1, i.e. measure all energies in units 
of half the band width. For finite T and not too large 
U the functional integral can be calculated numerically 
by Quantum-Monte-Carlo simulations Jl8|,^5| . Equ. (|J) is 
solved by iteration. All calculations were performed for 
U = 2 where the Neel temperature in zero field is close 
to its maximum value T^^ max = 0.10. 

The results for the magnetization m(H) and the stag- 
gered magnetization m st {H) |2(| are shown in Fig. 2. 
The metamagnetic behavior is clearly seen: at T = 
1/12 m(H) begins to show a typical "S-shape" which be- 
comes more pronounced at T = 1/16 (Fig. 2a). For these 
temperatures m s t is a continuous function of H (Fig. 2b), 
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which vanishes at the critical field H c as m s t ~ (H c —H) u . 
A least-square fit through five points (weighted with the 
error bars) near the transition yields v = 0.46 ± 0.05 for 
T = 1/12, which is compatible with a mean-field expo- 
nent v = 1/2. By contrast, at the lowest temperatures, 
T < 1/32 (Fig. 2c, d) the transitions in m and m s t are 
clearly discontinuous. Hysteresis is only weak and is not 
shown. The field dependence at intermediate tempera- 
tures 1/16 < T < 1/32 is more complex: m st is almost 
field-independent (i. e. m ~ 0) below some field H' c , de- 
creases sharply at H' c (such that m > 0), but vanishes 
only at a field H c > H' c in a continuous way. Although 
the error bars do not permit an unambiguous interpre- 
tation it seems that the order parameter decreases by 
two consecutive transitions: the first one, at H' c , being of 
1. order and the second one, at H c , of 2. order. 

The results for m(H) and m st (H) are used to construct 
the H — T phase diagram shown in Fig. 3. It displays all 
the features of Fig. lb. In particular, the 1. order line, 
defined by H' C (T) for T > Tce, continues into the ordered 
phase, separating two different AF phases: AF/9412107 
(where m ~ 0) and AF// (where m > 0). The position of 
its endpoint cannot, at present, be determined accurately 
(dotted line). We note that the ratio of Tqe and the Neel 
temperature T/v at H — is Tce/Tm — 0.3. 

The H — T phase diagram in Fig. 3 is strikingly similar 
to the experimental phase diagram of FeBr2 jl2|,[27j . In 
particular, the part of the 1. order line extending into 
the ordered phase clearly resembles the regime where 
de Azevedo et al. |l2| observed an anomalous, but non- 
critical, behavior in m(T) and x(T). It appears that 
there may well exist a truely critical feature in this part 
of the phase diagram. It would therefore be interest- 
ing to measure m s t (T, H) in FeBr2 by neutron scattering 
to see whether such a feature really exists. - The fact 
that the results shown in Figs. 2,3 are qualitatively very 
similar to the mean-field solution of the Ising model (|l|) 
with J' -C J [fipTl is surprising since at U = 2 (= band- 
width) the itinerant nature of the electrons cannot yet be 
ignored. On the other hand, the applicability of mean- 
field theory itself to FeBr 2 seems to be justified by the 
fact that in this material the AF superexchange involves 
20 equivalent sites in the two neighboring iron planes [p8| . 
We note that in FeBr 2 T CE /T N = 0.33 [|. 

In summary, we investigated the origin of metamag- 
netism in strongly anisotropic antiferromagnets start- 
ing from a microscopic model of strongly correlated- 
electrons, the Hubbard model, in a dynamical mean- 
field theory. The H — T phase diagram obtained for 
an intermediate on-site interaction U = 2 is qualita- 
tively very similar to that of FeBr 2 . Apart from cal- 
culations at U < 2, it will be interesting to calculate off 
half-filling. This will allow us to investigate the proper- 
ties of metallic metamagnets, such as the Uranium-based 
mixed-systems, for which a theory in terms of an itiner- 
ant electron model is mandatory. 
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FIGURE CAPTIONS 

1. Schematic H — T phase diagram for a) a typi- 
cal Ising-type metamagnet (TCP: tricritical point), 
b) the Ising model (Q) in mean-field theory with 
R < 3/5 (CE: critical endpoint, BCE: bicritical 
endpoint). Full lines: 1. order transition, broken 
lines: 2. order transition; AF: antiferromagnetic 
phase, P: paramagnetic phase. 

2. QMC-results, including error bars, for the magneti- 
zation m(H) and staggered magnetization m st (H) 
as obtained for the d — 00 Hubbard model with 
easy axis along H at n = 1 and U = 2 for two sets 
of temperatures. Curves are guides to the eye only. 

3. H— T phase diagram for the d = 00 Hubbard model 
with easy axis along H at n = 1 and U — 2 as 
constructed from the QMC-results for m(H) and 
m s t{H); same notation as in Fig. lb. 
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